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Abstract. Berry phases and quantum fidelities for interacting spins have attracted 
considerable attention, in particular in relation to entanglement properties of spin 
systems and quantum phase transitions. These efforts mainly focus either on spin 
pairs or the thermodynamic infinite spin limit, while studies of the multipartite case of 
a finite number of spins are rare. Here, we analyze Berry phases and quantum fidelities 
of the energetic ground state of a Lipkin-Meshkov-Glick (LMG) model consisting 
of three spin-^ particles (qubits). We find explicit expressions for the Berry phase 
and fidelity susceptibility of the full system as well as the mixed state Berry phase 
and partial-state fidelity susceptibility of its one- and two-qubit subsystems. We 
demonstrate a realization of a nontrivial magnetic monopole structure associated with 
local, coordinated rotations of the three-qubit system around the external magnetic 
field. 
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1. Introduction 

Recently, the study of various ground state properties of many-body systems has 
attracted a lot of attention in the context of quantum phase transition (QPT). QPTs, 
driven solely by quantum fluctuations, exhibit a dramatic change in the ground state 
at zero temperature under change of external parameters and is associated with a level 
crossing or avoided crossing between the ground state and first excited state |lj . A deeper 
understanding of QPTs has emerged from the fundamentals of quantum mechanics. In 
particular, when the external parameters of the Hamiltonian are varied, the response of 
the energetic ground state of the system has been analyzed in terms of Berry phase [2] 
and Bures-Uhlmann fidelities [3l H]. These quantities show nontrivial behavior related 
to points in parameter space where two or more energy levels become degenerate. 

For spin-like systems, energy crossings provide realizations of magnetic monopole 
structures; a fact that has opened up for simulations of magnetic monopoles in the 
laboratory [5]. These monopoles give rise to a magnetic flux that can be measured as 
a Berry phase being proportional to the area (solid angle) enclosed in the parameter 
space of the system. The monopole structure of interacting spins has been examined in 
several studies in the past P [3 El El [101 E] ■ 

The Berry phase for systems with several spin— | (qubits) has been addressed 
recently. Lei Xing [12] examined a three-qubit model with uniaxial qubit-qubit 
interaction and demonstrated that the corresponding Berry phase admits a solid angle 
interpretation, provided the couplings are added to the underlying parameter space. 
Williamson and Vedral fl3\ IH] found a nontrivial relation between the Berry phase 
of translationally symmetric multi-qubit states and their multi-partite entanglement 
properties. The behavior of the Berry phase in the thermodynamic limit has been 
studied for XY spin-chains [151 [161 [T71 [H] , the Dicke model [19], and the Lipkin-Meskov- 
Glick (LMG) system [20]. In this paper, we analyze the Berry phase for an analytically 
solvable, finite-size LMG type model [21] consisting of three spin— ^ particles. We 
examine the Berry phase structure of the energetic ground state of the LMG system 
and its subsystems using pure state [2] and mixed state [22] Berry phases. 

The fidelity [23] is an information theoretic measure that can be used to analyze the 
quantal properties of the ground state of spin systems. The utihty of this measure and 
the related fidelity susceptibility, has been explored in a number of studies [211 ISSl |26] , 
in particular in relation to QPTs. The fidelity susceptibility and the Berry phase are 
two complementary manifestations of the underlying geometry of the state space, as 
described by the quantum geometric tensor [271 [2H] • Recently, the concept of partial- 
state fidelity has been developed, which measures the fidelity of a subsystem along with 
the associated notion of partial-state fidelity susceptibihty [291 EO]- Here, we wish to 
examine the response of the three-qubit LMG ground state to parameter variation in 
terms of the fidelity susceptibility and partial-state fidelity susceptibility. 

The outline of the paper is as follows. In the next section, the three-qubit LMG 
model is described, the corresponding energy eigenvalues and eigenvectors are found. 
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and the ground states are identified. Sees. |3] and H] examine in detail the Berry phases 
and fidehty susceptibihties of the present LMG system and its subsystems. The paper 
ends with the conclusions. 



2. Three-qubit LMG model 



The LMG model of spin systems has found applications in Bose-Einstein condensates 
[31], statistical mechanics of mutually interacting spins [32], and entanglement theory 
[33| IM] . The LMG model describes a set of qubits (spin-|) mutually interacting 
through a XF-like term in the Hamiltonian and coupled to an external transverse 
magnetic field. The ferromagnetic version of the LMG Hamiltonian reads 



H{^,h-N) = --{Sl + ^Sl)-hS., 



(1) 



where 7 is an anisotropy parameter (7 = 1 corresponds to the isotropic LMG model), 
h is the strength of an external magnetic field in the z direction, and Sa = Ylk=i 
is the ath component of the total spin operator {h = 1 from now on) with o"^ = 
|0)(1| + |1)(0|,(T^ = -z|0)(l|+i|l)(0|, anda^ = |0)(0| - |1)(1| the Pauh operators of the 
kih qubit. In the present work, we examine the exactly solvable three-qubit (A^ = 3) 
case. 

The Hamiltonian for the three-qubit LMG system is given by 



1 



i/(7, h) = //(7, h- 3) = -- [alal + alal + alal + ^{ay^ + 0^ + aja^ 



h 



- [a] + al + af) , 



(2) 



where we have ignored the unimportant constant term — 1(1 +7)- In the computational 
basis {|000),|011),|101),|110),|111),|100),|010),|001)}, the Hamiltonian takes the 
block- diagonal form 



H{j,h) 



P{l,h) 






P(7,-/i) 



(3) 



where 



Ik 



-1(1-7) 



-1(1-7) 
-1(1-7) 

V 4(1-7) 



-1(1-7) 
-1(1 + 7) 



-1(1-7) \ 
-1(1 + 7) 
-1(1 + 7) 



(4) 



-1(1 + 7) -I 
-1(1 + 7) -1(1 + 7) 
and is the 4x4 null matrix. The block structure originate from the existence of 
the conserved quantity alafaf, which displays the fact that only states with same spin 
parity interact [35]. Due to the h ^ —h symmetry between the two P blocks of H, we 
may assume h > without loss of generality. To facilitate the diagonalization of H, we 
define energy functions S± and £, and mixing angle G according to 



f±(7,'>) 



7 



± 2^/9h^ + 3h{l + 7) + 1 - 7 + 72) 
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Soh,h)±AS{-f,h), 



tan 



v^(7 - 1) 



(5) 



6/i + 1 + 7 - 2 v/9/i2 + 3/^(1 + ^) + 1 _ ^ + y ■ 

Note, in particular, that tan [|6(7, /i)] (tan [i6(7, — /i)]) diverges (tends to zero) in 
the isotropic hmit 7 — ?■ 1. Thus, G(l, /i) = vr and —h) = 0. We may now write the 
eigenvalues E and orthonormalized eigenvectors \V) oi H in terms of £±.,£., and as 



E 



{+) 



eI-^ 



^+(7,/i): |y|+^) = -sin 



-9(7, /i) 



cos 



|000) 
1 000) + sin 



cos 



-6(7, /i) 



-9(7, /i) 



sm 




^(1011) -|110)), 
^(|011)-2|101)- 



|110)), 



1111) 



cos 



-9(7, -/i) 



-9(7, -/^) 



|111) + sin 



-9(7, -/i) 



= ^(|100)-|001)), 
= ^(|100)-2|010) + |001)), 

Oil) + |101) + |110)). 



\W), 

(6) 



where | W) = ® (Tx ® CTni:\ W) 

Alternatively, we may use that the total spin S"^ commutes with the LMG 
Hamiltonian [35], which implies that the eigensolutions may be labeled by the 
total spin. For instance, we may write the two types of ground states as 
|yi±)) = cos [i9(7,±/i)] If, ±1) + sin [i9(7,±/i)] 11,^1) with \S,M) being the 
common eigenvectors of S"^ and Sz- 



The lowest energy state is V^"^ or ' associated with energies E^_' and E. 



(+) 



n(-) 



Berry phase and fidelity susceptibility of the three-qubit Lipkin-Meshkov- Click ground stateh 



respectively, shown in FigHJ The crossing points he along the lines j ^ h^ = hci-y), 
which are determined by 

£_(7,M=^-(7,-M- (7) 
This yields the following two classes of solutions 

h^^ = 0, 7 arbitrary, (8) 

and 

= Ir (9) 

The energetic ground state in the low field strength regime, corresponding to h^c^ < h < 



he , is V_ ; in the high field strength regime, corresponding to /i > /ic , it is V_ . 

The energetic ground state may take any of the three main three-qubit forms: 
W, GHZ, and product states. Here, we identify the corresponding mixing angles and 
delineate the exact form of these ground states. The W and product forms are obtained 
for mixing angle being an integer multiple of tt, which may happen only in the isotropic 
case 7 = 1. Indeed, we found above that 9(1, /i) = tt and Q{l,—h) = 0, which implies 
the energetic ground state 

I^W) = /iW <h< hf\ 

= |iii), /i>/if (10) 

in the isotropic LMG model. The GHZ form Ui®U2® ^/g^dOOO) ± |111)), f/i, t/s, ^/g 
any one-qubit unitary operators, requires mixing angles that satisfy tan (|6) = ±\/3. 
This happens at 6 = 6(0,0) = ^ and 9 = 6(2, -|) = ^. Here, 6 = ^ corresponds 
to the two-fold degenerate ground state 

= a^(^a^® a,-^(|000) ± |111)), (11) 

i.e., GHZ form with Ui = U2 = = a^- Thus, the energetic ground state tends to a 
GHZ when approaching the origin in the (7, h) plane. The angle B = ^ yields a GHZ 
with Ui = U2 = U-i = ay. However, the corresponding states at (7, h) = (2, ±|) are 
yi^'*, neither of which being the energetic ground state. In other words, the ground 
state may be of GHZ form only at (7, h) = (0, 0). 



3. Berry phase 

3.1. Full system 

Here, we examine Berry phases [2j arising in adiabatic variation of the LMG 
Hamiltonian. For given 7 and h, let us consider the isospectral one-parameter 
Hamiltonian family 

H{-f, h- 0) = e~'^^'H{-f, h)e'*^% (12) 

where (p is slowly varying. Note that the unitary operator e~^'^^^, corresponding to 
coordinated rotation of the system around the z axis by an angle (p, preserves the 4x4 
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block-structure of H{-y,h). Thus, by preparing the system in the energetic ground 
state and by varying (p slowly, the system remains in the corresponding two-dimensional 
subspace. The state of the system may be represented by one of the double-valued 
eigenvectors 



-i(/.S',|y(±)^ 



(13) 



After completion of a 27r rotation around the z axis, corresponding to increasing from 
to 27r, we obtain the Berry phase in cyclic adiabatic evolution as [361 EZ] 

f = arg(yW(O)|yW(27r))+^/„''^(FW(0)|^yW(0))rf0 



(3, 



-2^(1 -cos [0(7, /i)]), h'i^'^ < h < h'^c^ 



(14) 



arg 



{Vl-^ (0) I K^-) (2^) ) + t C {V^~^ (0) I i;Vl-\cP)) d<l> 



h > h'c^ 



27r(l-cos[e(7,-/i)]), 

The absolute value |/3g| of the ground state Berry phase (3g is shown in Fig. [2l It should 
be noted that the Berry phase is defined modulus 27r, which implies that the 47r jump 
at the crossing point h^c^ = § the isotropic (7 = 1) LMG model that is visible in Fig. 
[2] cannot be detected experimentally. 

In order to understand the origin of the nontrivial two- level type Pg in Eq. (HM . 
we project the Hamiltonian -^(7, h; (p) onto two-dimensional subspaces spanned by 
{|000), \W)} for # <h< # and {|111), \W)} for h > hf\ Let = |000)(000| + 
and P*^^^ = |111)(111| + be the corresponding projection operators. 

-2|ooo)(W| + i|W)(ooo|, 



Furthermore, we define = |000)(1^| + |1^)(000|, 
and Si+^ = |000)(000| - as well as S^"^ = $ 



{+) 



ot 

0.0 



Figure 2. The absolute value \f3g\ of the ground state Berry phase (3g as a function of 
the isotropy parameter 7 and the magnetic field strength h. (This figure is in colour 
only in the electronic version) 



Berry phase and fidelity susceptibility of the three-qubit Lipkin-Meshkov- Click ground state! 



k = x,y, z. This yields the effective projected two-level ground state Hamiltonian 

' p(+)i7(7,/i,0)p(+) = E^{ci,h)P^^'^ 



+A£(7, h) sin [9(7, h)\ cos(2(^)S^: 



(+) 



Hg{-f,h,(j)) = < 



P^-^H{-i,h,(t))P^- 



+ sin [e(7,/i)]sin(2(/))S 

+ COS [e(7,/l)]Si+)'; , 



h^}^<h<hf\ 

£,(^,-h)P(-) + 



(+) 



(15) 



A£(7, -h) sin [6(7, -h)] cos(2(/))S^" 



-sin [6(7, -/i)] sin(20)S 



(-) 



+ COS [6(7,-/i)]S 



(-) 



h > h^c\ 

This describes a spin— ^ particle exposed to an effective magnetic field with strength 
A£^(7, ±/i) that takes the form 

= A£(7, ±/i) {sin [6(7, ±/i)] cos(20), ± sin [6(7, ±h)\ sin(2(/)), cos [6(7, ±h)]} 



= A^(7,±/i)n(7,±/i,0). 

Here, we have introduced the unit vectors 11(7, ±/i, (p) that rotate with twice the spin 
rotation angle and make polar angles 6(7, ±/i) with the effective z axis. The sign 
difference in the expression for the two types of ground state Berry phases in Eq. f|T4l) 
originates from that B^^^ rotate in opposite direction. It is visible that the origin 
of the nontrivial ground state Berry phase is a monopole sitting at the point where 
A£^ = 0. This happens at (7, /i) = (1, — |) in the low field regime {h^P < h < h^iP) 
and at (7, /i) = (1, |) in the high field regime {h > h^c^)] for 7 7^ 1 there is an avoided 
crossing at /i = — |(l + 7) (/i = |(l + 7)) corresponding to a minimal energy difference 
2A£^ = -igll — 7I in the low (high) field regime. The Berry effective gauge field takes 
the magnetic monopole form 

1 n(7, ±/i, (;/)) 



(16) 



^eff 



(17) 



2[A£:(7,±/i)]' 

and the Berry phase shown in Fig. |2] is the flux of -Bj^g through any surface enclosed 
by the curve traversed in parameter space (A£^, 6, 2^), where ^£ and 6 are determined 
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(2) 

by 7 and h. We may therefore interpret the jump across the crossing hne 7 /ic as 
an interplay between a jump in polar angle G and that the two types of ground state 
feel monopoles sitting at different points in the (7, h) plane. 



3.2. Subsystems 



An interferometer experiment to detect the Berry phase could be set up for one or two 
of the qubits. As the states of the subsystems in general are mixed, the corresponding 
Berry phases would coincide with the mixed state geometric phase in Ref . |22] , applied 
to adiabatic evolution. Here, we examine the behavior of these mixed state Berry phases 
in the LMG system. 

We calculate the subsystem Berry phases under slow rotation around the z axis. 
To this end, we need the reduced ground states p^"^^ and of the one- and two-qubit 
subsystem, respectively. Taking into account the translational symmetry of the ground 
states vd^'', these marginal states for any of the qubit or qubit pair read 

pW = l[i + r(7,±/iK] 



(18) 



and 



(±) = - 



[l + r(7,±/i)]|^r)W 



+ -[l-r(7,±/i)]|^2 



(19) 



respectively. Here, 
r(7, /i) 



[1 + 2 cos 6(7, /i)] 
1 

'2 + cos [6(7, /i)] 
1 



1 cos 



^f2 

1 
72 



h COS [6(7,-/1)] 
|01) + |10)' 



20(7, /i) 



1 00) + sin 



111) 



1 11) + sin 



5i+i®<Tz)| / (±)\ 



li^l^') of f?^^''(0) corresponding to its 



which define the eigenvectors e~^2'^ 
nonzero eigenvalues. 

We notice that the one-qubit Berry phases vanish since the corresponding 
reduced density operators are diagonal in the |0), |1) basis and thereby commute 



with e 



On the other hand, the two-qubit Berry phases may be non-vanishing. 



To see this, we note that the reduced two-qubit density operator reads ^'•^^(0) = 
g-i|0{a,®i+i®a,)^(±)gii0(<x,®i+i®<T,)_ ^-^^^ ^(^^0) ^ g^^\ which opens up for 

nontrivial mixed state Berry phase F*^^^ of the two-qubit states. 
By applying [22] to a cyclic adiabatic evolution, we obtain 



r(±) 



arg 



(21) 
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Figure 3. The absolute value \Tg\ of the two-qubit Berry phase Tg of the reduced 
ground state as a function of the isotropy parameter 7 and the magnetic field strength 
h. (This figure is in colour only in the electronic version) 



where p^it^ and /Sj^^^ are the density operator's eigenvalues and eigenstate Berry phases, 
respectively. For cyclic evolution, we obtain from Eq. (12T]) the two-qubit geometric 
phase Tg of the ground state as 

f rW = arg|(2 + cos[e(7,/i)])e*'^i'" + (1 -cos[e(7,/i)])e^'52'''} 



arctan 



(2+cos[e(7,fe)])sin/3^+' 



l-cos[0(7,h)] + (2+cos[e{7,/i)])cos /^J^ 



h^c^ <h<h. 



{+) ' "-c 



(2) 



(22) 



r(-) = arg I (2 + cos [6(7, -h)]) e^^^i"' + (1 - cos [6(7, -h)]) e*'^^ ' | 



arctan 



(2+cos[e(7 -fe)])sm/3^ ' 



where we have used that j3'^'^ = since are eigenvectors of cr^ ® 1 + 1 ® cr^. Here, 



l-cos[e(7,-/i)]+(2+cos[e(7,-/i)]) cos 01 

(±) 



(-) ' 



^(.) ^^2.,^ -cos [9(7, ±/.)] 



(23) 



2 + cos [6(7, ±/i)]' 

which is /3^^^ quenched by a factor {2 + cos [6(7, ±/i)]}~^. The absolute value IF^I of 
the two-qubit Berry phase Tg of the ground state is shown in Fig. [3l 

Note that the relative phase arg Tr (^e"*^'^('^^'^^+^®'^^)^(=^) j [22] between the 
initial and final two-qubit states is measurable for all (7, /i), as the visibility 



cannot vanish for these states. On the other hand, the 
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geometric part Tg of this relative phase is not always well-defined; in fact, it is undefined 
if ^7*^^) has nonzero degenerate eigenvalues [2^, i.e., when cos0 = — |. Thus, the mixed 
state Berry phase is undefined precisely when the three-qubit ground state is of GHZ 
form, i.e., at (7, /i) = (0,0), see Eq. (fTTI) . Note also the oscillatory behavior of Tg 
that is visibility in Fig. |3] when approaching the degeneracy. These oscillations may 
be understood from the fact that the mixed state Berry phase in unitary evolution is 
known to vary more rapidly in the vicinity of a degeneracy point of the corresponding 
density operator [3H]. Finally, just as the three-qubit Berry phase factor, the two-qubit 
Berry phase factor is smooth across the crossing point at /i = | in the isotropic LMG 
model. Indeed, vanishes for all h as is an integer multiple of 27r when 7 = 1. 



4. Fidelity 

4.1. Fidelity susceptibility 

Fidelity is a measure of similarity between different quantum states and is therefore 
expected to be sensitive to abrupt changes in the ground state properties in many-body 
systems. This has triggered work to use fidelity measures in the context of quantum 
critical phenomena |23l |271 |28l [211 ISSl |26]. Here, we analyze the fidelity susceptibility 
[23] in the present three-qubit LMG system. 

An abrupt change in the LMG model system can be induced by slowly tuning the 
external magnetic field h across the crossing value h^^ = |a/7 at fixed 7. The fidelity 
susceptibility Xgil^h) is taken to measure the response of the ground state to such 
variations in h. A convenient form for this Xgili h) may be found by writing the LMG 
Hamiltonian as 

H{-f,h) = Ho + hHi, (24) 

where Hq = —^(^Sl + ^ySy) is independent of h and Hj = —S^ is the driving 
Hamiltonian, the relative strength of Hq and Hj being controlled by h. Let \Vg{'-f,h)) = 
iV^i^"*) be the normalized ground state of H{'y,h) and Eg{'y,h) the corresponding 
ground state energy. The fidelity susceptibility Xgili h) of Vg is defined as the leading 
nontrivial contribution in dh to the fidelity Fg between the ground states |K;(7; h)) and 
\Vg{'y, h + Sh)), according to 

F,(7, h, Sh) = I {Vg{^, h)\Vg{^, h + 6h))\ 

1 

r 

where we may note that Xg is independent of the arbitrary parameter dh. By expanding 
to second order in 6h and using the form of i^(7, h), we obtain 

J (26) 

where |T4(7, h)) and i?n(7, h) are eigenvectors and eigenvalues, respectively, of H{j, h). 



-X';{l,h)6h' + ..., (25) 
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By inserting Eq. into Eq. (B^ and using = —S^, we obtain 



\{vi-'>\{-s.)\v[-^)\ 



sin^[e(7,fc)] 
sin^[e(7,-fe)] 



(27) 

h > hf\ 



We may note that 5"^ does not couple the degenerate states V^'^'^ and ¥2"^^^ to any of 
the two candidate ground states, since ^ and belong to different S'^ eigenvalues. 
The fidelity susceptibility Xg(7, h) is shown in Fig. HI 




Figure 4. Fidelity susceptibility Xg of ground state as a function of the isotropy 
parameter 7 and the magnetic field strength h. (This figure is in colour only in the 
electronic version) 

The fidelity susceptibility of the LMG ground state vanishes identically for 7=1 
as G(l,/i) = or TT. therefore shares the behavior of the Berry phase factors e*^^ 
and e*'"^ in that it is smooth across the crossing point in the isotropic (7 = 1) case. 
Furthermore, x^ili ^c^^ ) decreases monotonically towards zero as a function of 7 when 
7 increases. This can be explained by noting that sin^(G) decreases monotonically when 
9 increases from 9(0,0) = ^ to 9(l,/ii^^~) = vr. On the other hand, X^ilihf^^) has 
a local maximum since to decrease 9(0,0) = ^ to 9(1, /ii^^^) = vr one must pass the 
intermediate angle 9 = ^ at which sin 9 has its maximum. Furthermore, we note 
that the fidelity susceptibility is singular close to the degeneracies at (7, /i) = (1,±|) 
which corresponds to the locations of the effective magnetic monopoles and where the 
adiabatic approximation breaks down. The singular behavior expresses the fact that 
small variations in the parameters may cause transitions between the two orthogonal 
states that cross at these points. 
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4-2. Partial-state fidelity susceptibility 

Partial-state fidelity susceptibility has been developed to deal with the response of a 
subsystem s to the driving Hamiltonian [29| 130] . It is defined as the leading nontrivial 
contribution of the Bures-Uhlmann fidelity Fg.g [3l H] of two marginal ground states 
Pg{^, h) = Trp|y,(7, h))iyg{^, h) I and p,(7, h + 5h) = IVp|V'g(7, h + (^3(7, h + 5h)l 
Tip being partial trace over one or two of the qubits. Explicitly, 

F,.g(7, /i, 6h) = sj pg{-i,'h)pg{-f,h + 6h)^J pgin.h) 



= l--Xs;,{l,h)6h' + ... 
defines the partial-state susceptibility Xs n of ^he energetic ground state [29| 130]. 



(28) 




1.0^00 



Figure 5. One-qubit partial-state fidelity susceptibility Xi-g of the LMG ground state 
as a function of the isotropy parameter 7 and the magnetic field strength h. (This 
figure is in colour only in the electronic version) 



Let US first consider the one-qubit partial-state fidelity Xi-g with respect to 
variations of h. The marginal ground state p'-^^ of any of the three qubits is diagonal in 
the fixed |0), |1) basis. This implies that only changes in the purity parameter r{j,h) 
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contribute to Xi-n- explicit calculation yields 



2 l-[r(7,h)]2 idh^i'^^ ^)] 

r^er-Y h)V 



4-2 cos ( 



hi^^ <h< hf\ 



Xi{l,-h) = lTZT^f-Kw[iir{-f,-h)] 



(29) 



sin^ [e(7,-fe)] 



4-2cos[0(7,-h)]-2cos2[e(7,-/i)] 



h > hf\ 



The one-qubit partial-state fidelity susceptibility Xq(7, h) is shown in Fig. [51 




0.5 h 



Figure 6. Two-qubit partial-state fidelity susceptibility X2-g of the LMG ground state 
as a function of the isotropy parameter 7 and the magnetic field strength h. (This 
figure is in colour only in the electronic version) 

The two-qubit partial-state fidelity X2;g originates from changes in the purity of 
Q^'^^ and in the parameter dependent eigenvector l'?/'!^^) of q'^'^\ The relevant purity 
parameter r(7, /i), which implies that the contribution to X2 g from the purity coincides 
with the one-qubit partial-state fidelity susceptibility Xi-g- The additional contribution 
related to the change of \tp^^^) equals the corresponding pure state fidelity susceptibility, 



Berry phase and fidelity susceptibility of the three-qubit Lipkin-Meshkov- Click ground statelA 



weighted by the probabihty Ml + ^^(7, Exphcitly, we have 



X2h,-h) 



2 l-[r(7,/i)]^ 



hl^^ <h<h 



(2) 
c ) 



1 1 r^, 

2 l-[r(7,-fc)]2 Vdh' 



(30) 



+ [1 + ^7,-/^)] 



V 

The two-qubit partial-state fidehty susceptibihty X2;g(7i ^) is shown in Fig. [61 

Both the one- and two-qubit partial-state fidelity susceptibilities behave similarly 
as that of the full ground state Vg-. both Xi-gili h) and X2 g(7i ^) vanish in the isotropic 
(7 = 1) case and there is a similar dependence on 7 close to the crossing line 7 1— )■ /ic . 
Furthermore, by comparing Figs. HI El and El we notice that the pure state fidelity 
susceptibilities is typically larger than the partial-state fidelity susceptibilities. It is 
apparent that X2;g(7i ^) — Xi-gili h) since the second term in the right-hand side of Eq. 
fl5U]) . This may be interpreted to be a consequence of the loss of purity for each qubit 
that is traced out. 



5. Conclusions 



A detailed characterization of the ground state of a three-qubit Lipkin-Meshkov-Glick 
(LMG) type model has been given. We have calculated Berry phases for the three-qubit 
state as well as for the reduced two-qubit state in the case of local, coordinated 2% 
rotation around the axis of the external magnetic field. We have identified an underlying 
two-level structure of the three-qubit Berry phase and found the relevant magnetic 
monopole distribution. The energetic ground state of the model is of GHZ-type if the 
external field and the isotropy parameter both vanish. The reduced two-qubit state at 
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this point in parameter space is two-fold degenerate and separable, from which follows 
that the corresponding mixed state Berry phase is undefined. The three- and two-qubit 
Berry phases vanish modulus 27r in the isotropic LMG model. 

We have calculated the fidelity susceptibility and the one- and two-qubit partial- 
state fidelity susceptibility for the LMG model. These fidelity susceptibilities all behave 
similarly, but decreases in size for each qubit being traced out. We have found that 
the fidelity susceptibilities all vanish in the isotropic LMG model. Analogously to 
the fidelity susceptibility and Berry phase in the pure state case, the partial-state 
fidelity susceptibility and the Uhlmann holonomy [39] measure the geometry of the 
space of mixed quantum states. This observation makes it natural to ask whether 
the Uhlmann holonomy may yield further insights into the ground state properties of 
interacting spin models. Paunkovic and Rocha Vieira [40] have found a rich structure 
in the Uhlmann holonomy for thermal states in the Stoner-Hubbard and BCS models. 
A similar calculation of the partial-state holonomy seems pertinent in relation to the 
present work. 

We hope that the analysis presented in this work may trigger investigations of 
few-qubit models to explore further their effective magnetic monopole structure and its 
associated state space geometry. 
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